Density Matrix Renormalization Group study of ^^Cr and ^^Ni 
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We discuss the development of an angular-momentum-conserving variant of the Density Matrix 
Renormalization Group (DMRG) method for use in large-scale shell-model calculations of atomic 
nuclei and report a first application of the method to the ground state of ^^Ni and improved results 
for ^^Cr. In both cases, we see a high level of agreement with the exact results. A comparison of 
the two shows a dramatic reduction in the fraction of the space required to achieve accuracy as the 
size of the problem grows. 
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One of the foremost challenges confronting nuclear 
physics today is the systematic study of medium-mass 
and heavy nuclei using the shell model. Even with the 
drastic truncation achieved by limiting the active parti- 
cles to one or at most a couple of valence shells outside a 
doubly magic core, the size of the resulting space still ex- 
ceeds storage capabilities of the currently available com- 
putational resources for all but fairly light nuclei. This 
opens up the need for innovative truncation strategies. 
The Density Matrix Renormalization Group (DMRG) 
method has had outstanding success dealing with low- 
dimensional quantum lattice problems [ij. The method 
was later extended to finite fermi systems, where it has 
been applied with impressive success to the description 
of small metallic grains Qito two-dimensional electrons 
in strong magnetic fields [3[, and to problems in quan- 
tum chemistry \A] . This suggests that it might also prove 
useful in the description of another finite fermi system, 
the atomic nucleus. In this work we present our first re- 
sults for ^^Ni, the most ambitious test we have considered 
to date, and confirm its usefulness for nuclear structure 
calculations. 
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The usual DMRG algorithm begins by partitioning the 
complete Hilbert space for a given problem onto a set of 
lattice sites, with each site admitting a set of basis states. 
The untruncated problem is recovered by considering the 
complete set of states within all sites. Since the full space 
is typically too large for exact treatment, the DMRG 
method treats the problem iteratively, by successively 
adding sites to those already treated and implementing a 
truncation based on density matrix considerations. This 
is illustrated schematically if Figure 1. 

Assume that there are n sites on the lattice and that r 
sites to the left have already been treated, defining what 
we call the system block B^- Assume further that within 
this block we have m states and the matrix elements of 
all suboperators of the hamiltonian for these states. The 
next step would be to enlarge by adding the r + 1** 
site, producing the enlarged block Br+i = Br®Lr+i- We 
then couple this to a medium, which involves information 



FIG. 1: Schematic illustration of the DMRG growth proce- 
dure. Shown are all the sites in the lattice chain, those of the 
system block, those added in the enlargement process, those 
of the medium, and their organization in the superblock. 



on all of the remaining sites, producing the Superblock 
Br ® Lr+i ® M, representing the entire system. The 
basic idea of the DMRG method is to truncate the en- 
larged block to its optimum m states, namely those that 
are the most significant contributors to the Superblock 
ground state, and to renormalizc all operators to act in 
this truncated space. 

We will assume for now a product space description 
and denote the states of the enlarged block as \I)b,.+i — 
|*)s^|j)r--i-i, where i span the states of the system block 
and j those of the added site. The Superblock ground 
state can then be written as 
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where k spans the states of the medium (M). 

The reduced density matrix for the enlarged block in 
the Superblock ground state is obtained by contracting 
over the states of the medium, namely 
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If we diagonalize this reduced density matrix and main- 
tain the eigenstates associated with its m largest eigen- 
values we are guaranteed to have found the m most 
important components of the enlarged block in the Su- 
perblock ground state. 
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The basic idea then is to systematically grow the sys- 
tem block by adding lattice sites and then at each stage 
to truncate to the m most important states obtained in 
this way. At each stage we must, as just noted, transform 
all hamiltonian suboperators to the m— dimensional trun- 
cated space, as this provides required information for its 
subsequent enlargement. The enlargement process just 
described is repeated over and over, sweeping in both di- 
rections through the set of lattice sites until convergence 
in the ground-state energy is obtained. The calculations 
can then be redone as a function of m until acceptably 
small changes with increasing values are obtained. 

The usual DMRG procedure, as just outlined, works 
with product states and is thus equivalent to an m- 
scheme approach in the nuclear context. This has the 
drawback that conserving angular momentum symmetry 
becomes difficult. To avoid this problem in applications 
to nuclei, Dukelsky and Pittel [5| proposed the use of 
an angular-momcntum-preserving variant of the DMRG, 
called the JDMRG. This method, which is an example of 
a non-Abelian DMRG algorithm was first applied to 
nuclei in the context of the Gamow Shell Model Q and 
then subsequently applied to the traditional nuclear shell 
model by Pittel and Sandulescu @ . 

In the JDMRG, we work in a coupled (or J-scheme) 
basis. As a consequence, we must now calculate reduced 
matrix elements of the various suboperators of the Hamil- 
tonian, namely 

al{a\a])\{a%)\[{a\a]Y~ak]\ 
[{a\a])\hk~ai)^]° + h.c. 

An important feature of nuclei is that they contain 
two types of particles, neutrons and protons. This leads 
to the question of how the associated orbitals should be 
arranged on the lattice, intertwined or at opposite ends. 
We have found it useful to maintain a separation between 
the neutron and proton blocks, leading to what we call 
a three-block JDMRG algorithm, for reasons to be made 
clear a bit later. 

The nuclei for which we will present results here are 
^*Cr and ^^Ni. In both we assume a doubly magic '*°Ca 
core with the remaining nucleons distributed over the 
orbitals of the 2p — If shell. We now spell out in a bit 
more detail the steps of the JDMRG procedure that were 
implemented in our treatment of these two nuclei. 

Setting the order of orbits: 

As just noted, we have chosen to define the lattice chain 
such that the neutron and the proton orbits lie on its 
opposite ends. This still leaves several possible options, 
however. We have found the optimal order to be the one 
shown in Figure 2, where we use the notation nj and pj to 
denote a neutron orbital with angular momentum j and 
a proton orbital with angular momentum j, respectively. 
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FIG. 2: The order of orbits used in the calculations described 
in the text. 
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FIG. 3; Two consecutive enlargement steps during the 
warmup phase of the DMRG calculations described in the 
text. In the first, the third neutron orbit from the left is 
added; in the next, the third neutron orbit from the right is 
subsequently added. 

This order places maximally entangled orbits next to 
one another, as is found to be important from quantum 
information considerations Q . 

Warmup Phase(The first enlargement stage): 

This phase of the DMRG procedure is aimed at defining 
a good initial set of states to use in each block. In our 
calculations, we follow the procedure schematically illus- 
trated in Figure 3, whereby we begin by adding two orbits 
at one end of the chain of sites, followed by a subsequent 
enlargement at the other end, followed by gradual further 
enlargements from the two ends of the chain. In this ini- 
tialization phase, the last truncated block on the opposite 
end of the chain is used as the medium. As should be 
clear, the Superblock so obtained does not contain all 
of the orbits. Hence the Superblock must be diagonal- 
izcd for all relevant neutron and proton numbers in the 
J = 0+ space, to permit some particles to reside in or- 
bits not included. A mixed density matrix is constructed 
from the corresponding eigenstates which is then used 
to truncate the system to the m most important states 
for each block of orbits. This process is repeated until 
the middle of the chain, at which point all neutron and 
proton orbits have been treated. 

At this point we have an initial choice for the m most 
important states for each relevant block. 

The Sweep Phase( Successive enlargement stages): 

We now wish to iteratively update information on each 
block, by gradually sweeping through the orbits and us- 
ing information initially from the warmup stage and later 
from the previous sweep stage to define the medium. We 
start the sweep phase by adding the orbits from the cen- 
ter of the chain towards the end. The medium now con- 
sists of one of the like-particle blocks (Mi) - constructed 
either during the warmup or the previous sweep, and the 
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FIG. 4: Three consecutive steps in the enlargement of the 
neutron block from the left during the sweep phase. 



unlike-particle block (M2). The fact that our medium 
contains two components is why we refer to this as the 
three-block algorithm. The complete growth of the neu- 
tron and proton blocks is repeated in both directions until 
convergence is attained. Several consecutive steps in the 
neutron block enlargement during a given sweep stage 
are illustrated schematically in Figure 4. Because of the 
great care taken in the warmup initialization, only three 
to four sweeps are typically needed for the process to 
converge. 

Convergence with m: 

The above steps are repeated with increasing m until 
satisfactory convergence in the ground-state energy is 
obtained. 

We present here test results for ''^Cr and ^^Ni, using 
the three block algorithm described above. We have used 
the KB3 matrix elements [13] used in Refs. 11, 13], with 
the same set of single-particle energies used there as well. 
Excited states for ^^Ni were not included in the analysis 
since they are not available for the hamiltonian that we 
used. 



TABLE I: Calculated ground-state energies in MeV as a func- 
tion of m for "^^Cr. The maximum dimension encountered in 
the sweep process is also given. 



m 


Egs 


Max Dim 


40 


-32.698 


1,985 


60 


-32.763 


2,859 


80 


-32.788 


3,765 


100 


-32.817 


4,494 


120 


-32.840 


6,367 


140 


-32.890 


8,217 


160 


-32.902 


9,978 


180 


-32.944 


11,062 


200 


-32.947 


12,076 


Exact 


-32.953 


41,355 



The results for the ground-state energy of ''^Cr are 
shown in Table L The full shell-model space of ^*Cr con- 
sists of 1,963,461 states of which 41,355 have J"" = 0+. 

The exact result for the ground state is -32.953 MeV. 
Though the JDMRG results converge monotonically with 



TABLE II: Calculated energies in MeV for low-lying excited 
states in ''^Cr, with dimensions shown in brackets. 



m 






E.+ 


E,,+ 




( Dim') 


( Dim\ 


( Dim\ 


( Dirn\ 


1 nn 




'\n Qfi 

ou. yu 


-9Q 1 fi 


97 Q7 




003"! 


too 5 Z(U -L ) 


('38 fi59') 


(44 Q4"l 


120 


-32.01 


-30.93 


-29.20 


-28.06 




(OS p,77'\ 


(AO ')'\A\ 


(AK r\<\A\ 


(a 'ia'7'\ 


140 


-32.04 


-30.98 


-29.26 


-28.15 




fSfi 7np^ 




1 oz 5 ufjyj 1 


rs 91 7~l 


1 fin 




-0 i. .KJ'-i: 




-28.29 




(AA MR) 


(63 4871 


f63 537") 


(9 9781 


180 


-32.13 


-31.09 


-29.43 
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(50,030) 


(72,616) 


(74,346) 


(11,062) 


200 


-32.13 


-31.10 


-29.47 


-28.48 




(54,891) 


(81,249) 


(85,168) 


(12,076) 


Exact 


-32.15 


-31.13 


-29.55 


-28.56 




(182,421) 


(246,979) 


(226,259) 


(41,355) 



TO to this value, we need a substantial portion (about 
25%) of the full space to achieve agreement to within a 
few keV of the exact energy. 

It should be noted that the results reported in Ref. 
[§] were for a different order of single-particle levels, 
explaining why the ground-state energies and maximum 
dimensions listed in that paper are different from those 
reported here for the same values of m. 

The results for low-lying excited states are provided 
in Table 2. To obtain these results, we construct the 
relevant hamiltonian matrices using the optimum block 
structures obtained at the point of the lowest ground 
state energy. Note that we do not target the excited 
states themselves in constructing the reduced density ma- 
trix used for the truncation. As a result, the convergence 
is not as rapid as for the ground state, but nevertheless 
quite acceptable. 



TABLE III: Calculated ground-state energies in MeV as a 
function of m for ^^Ni. The maximum dimensions encoun- 
tered in the sweep process are also given. 



m 


Egs 


Max Dim 


80 


-78.351 


72,023 


100 


-78.363 


83,773 


120 


-78.372 


102,690 


140 


-78.376 


119,797 


160 


-78.390 


136,073 


180 


-78.393 


162,019 


200 


-78.399 


192,878 


Exact 


-78.46 


15,443,684 



The results for ^^Ni ground state are tabulated in Table 
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3. The size of the exact space in an angular momentum 
basis is 15,443,684. We have carried out these calcula- 
tions with the same order of single-particle orbits as for 
**Cr. The fraction of space required for meaningful con- 
vergence is dramatically reduced for ^^Ni compared to 
*^Cr. With about 1% of the space we are able to get 
within around 60 keV of the exact results [l^. 

Earlier calculations for ^^Ni using the DMRG by Pa- 
penbrock and Dean [3], although done in the m-scheme, 
were unable to obtain the ground state energy to better 
than 400 keV. 

An interesting feature of our results for both *^Cr and 
^^Ni is the absence of an exponential falloff in the con- 
verged ground state energy as a function of the number 
of states retained. Since exponential behavior is often 
used to extrapolate to the actual ground state energy, it 
is worth commenting on why it does not occur here. 

A key difference between our JDMRG algorithm and 
the usual DMRG is that our lattice sites vary significantly 
in size, as they represent complete single-particle orbits. 
Often, when an orbit is added, no truncation is required 
at that step. When the increment in m reaches a large 
enough value that truncation is avoided, we find a sharp 
drop in the ground state energy for that m value. For 
lower values, truncation always arises at that step in the 
sweep and the falloff is less severe. As a result, the energy 
tends to fall off in spurts as a function of m rather than 
as a smooth exponential, with larger orbits producing 
a more jagged behavior. If, however, we carry out our 
calculations for a wide enough range of m values, a fit 
can still be carried out in the presence of such behavior. 
Later, we discuss another way of enhancing exponential 
falloff in the JDMRG. 

In this work, we have described our efforts aimed at 
developing the DMRG method as a dynamical truncation 
strategy for large scale shell model calculations of atomic 
nuclei. Following a brief description of the usual DMRG 
procedure, we discussed a specific three-block angular- 
momentum-conserving algorithm for carrying out such 
calculations. We presented results for two nuclei '^^Cr 
and ^®Ni. Although a large fraction of the space was 
required to achieve results of high accuracy for "^^Cr, the 
fraction of the space required for ^^Ni was much smaller. 
This bodes well for the future usefulness of the method 
for even larger problems. 

There are several issues that we plan to explore in the 
near future. On the one hand, we would like to extend 
our test analysis to a broader range of nuclei, to more 
meaningfully explore the requisite fraction of the space 
required for convergence. We also plan to redo our calcu- 
lations for ^"^Ni to use the GXPFlA interaction [li], for 
which exact results are available not only for the ground 
state but for the states in the first deformed band as 
well [l^l- This will help us to better assess the ability 
of the method to accurately treat excited states as well, 
even when only the ground state is being targeted in the 



iterative truncation process. 

Finally, to be able to study even heavier nuclei, it will 
be necessary to include orbits with even larger angular 
momenta. For example, in nuclei beyond ^®Ni the I59/2 
orbit becomes increasingly more important. The larger 
the orbit the greater is the computational strain it im- 
poses on the iterative growth process of the DMRG. A 
possible way to get around this is to split such orbits 
into two or more parts, while still maintaining exact an- 
gular momentum conservation throughout. We plan to 
test some possible approaches for splitting such orbits in 
the context of the smaller I/7/2 orbital and to compare 
the results that emerge with those we have already ob- 
tained. Hopefully this will enable us to then turn to even 
larger orbits and to continue our applications up in nu- 
clear mass. Use of smaller lattice sites may also facilitate 
a more rapid transition to exponential behavior in the 
ground state energy. 
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